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Abstract 

The influence of quantized electromagnetic fields on a nonrelativistic charged particle moving 
near a conducting plate is studied. We give a field-theoretic derivation of the nonlinear, non- 
Markovian Langevin equation of the particle by the method of Feynman- Vernon influence func- 
tional. This stochastic approach incorporates not only the stochastic noise manifested from elec- 
tromagnetic vacuum fluctuations, but also dissipation backreaction on a charge in the form of the 
retarded Lorentz forces. Since the imposition of the boundary is expected to anisotropically modify 
the effects of the fields on the evolution of the particle, we consider the motion of a charge under- 
going small-amplitude oscillations in the direction either parallel or normal to the plane boundary. 
Under the dipole approximation for nonrelativistic motion, velocity fluctuations of the charge are 
found to grow linearly with time in the early stage of the evolution at the rather different rate, 
revealing strong anisotropic behavior. They are then asymptotically saturated as a result of the 
fluctuation-dissipation relation, and the same saturated value is found for the motion in both 
directions. The observational consequences are discussed. 
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I. INTRODUCTION 



The manifestation of vacuum fluctuations can be visualized through a mechanical effect 
on a materialized body. One of the most celebrated examples is the attractive Casimir force 
between two parallel conducting plates [ij • This phenomenon in general can be characterized 
by fluctuations under geometric constraints such that its spectrum on long wavelength modes 
is modified. Nevertheless, this induced-force effect can also be probed through the coupling 
to a test particle. For example, consider an atom in its ground state as a test particle 
located near a perfectly conducting plate. The atom then experiences a position-dependent 
energy shift due to the boundary effect on vacuum fluctuations from which to give rise to 
an attractive Casimir- Polder force toward the plate Q]. Thus, the presence of the boundary 
is expected to anisotropically change vacuum fluctuations. In this paper a charged test 
particle serves as a probe to understand the nature of electromagnetic vacuum fluctuations 
by observing their effects on the test-particle's trajectory. 

When a charged particle interacts with quantized electromagnetic fields, not only do the 
expectation values of fields determine the mean trajectory of the charge, but the accom- 
panying quantum fluctuations also drive the charge into a zig-zag motion. The dynamics 
of the particle and field interaction has been studied quantum-mechanically in the system- 
plus-environment approach [3]. We treat the particle as the system of interest, and the 
degrees of freedom of fields as the environment. The influence of fields on the particle can 
be obtained with the method of Feynman- Vernon influence functional, by integrating out 
field variables within the context of the closed-time-path formalism |4|, |5|. The Langevin 
equation can then be derived by ignoring the intrinsic quantum uncertainty of the particle, 
which is assumed to be much smaller than the resolution of the position measurements. This 
stochastic approach incorporates both dissipative backreaction arising from the interaction 
with fields, and a stochastic noise owing to the quantum fluctuations of fields. In particular, 
the non- uniformity of the charge's motion will result in radiation that backreacts on itself 
through the electromagnetic self-force. The stochastic noise, which encodes the influence of 
quantum statistics of fields, drives the charge into a fluctuating motion p, |7|. Furthermore, 
the noise-averaged result reduces to the known Abraham-Lorentz-Dirac equation with the 
self-force given by a third-order time derivative of the position as expected {(J. 

The anisotropy of electromagnetic vacuum fluctuations in the presence of the conducting 
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plate has been studied via an interference experiment of the electrons, and is manifested in 



the form of the amp 
the previous article 



itude change and phase shift of the interference fringes 



In 



111 ], we employ the method of influence functional, and obtain the evo- 
lution of the reduced density matrix of the electron with self-consistent backreactions from 
quantized electromagnetic fields. Under the classical approximation with prescribed elec- 
tron's trajectories, it is shown that the modulus of the exponent in the influence functional 
describes the change of the interference contrast, and its phase results in an overall shift of 
the interference pattern. It is also found that the presence of the boundary anisotropically 
modifies the contrast of the interference fringes. In Ref. [3], the Brownian motion of a 
charged particle coupled to electromagnetic vacuum fluctuations near a perfectly conduct- 
ing plate is studied for the case that the particle barely moves; thus its dissipation effects is 
ignored. The behaviors of velocity fluctuations are shown different for the particle's motion 
in the directions perpendicular and parallel to the boundary plane. In this paper, we wish 
to further explore the anisotropic nature of vacuum fluctuations due to the boundary by 
the nontrivial motion of the charged particle where dissipative backreaction is incorporated 
in a consistent manner. In the presence of the boundary, we expect that radiation emitted 
by the charge in nonuniform motion should be bounced back and then impinge upon the 
original charge at later times. This will give rise to an additional retardation effect, which 
in turn results in a non-Markovian evolution of the particle. Thus, its mean trajectory will 
be altered. Here we will apply red the approach of influence functional from which to derive 
the Langevin equation beyond the mean-field approximation. The trajectory fluctuations off 
the mean value driven by the stochastic noise will be studied with dissipation backreaction 
taken into account in a way that a underlying fluctuation-dissipation relation is obeyed. 

Our presentation is organized as follows. In Sec. [Ill we introduce the closed-time-path for- 
malism to describe the evolution of the density matrix of a nonrelativistic charge coupled to 
quantized electromagnetic fields. We trace out field variables to obtain the coarse-grained ef- 
fective action from which the Langevin equation is derived. The resulting stochastic Lorentz 
force, which can be cast into a gauge invariant expression will be discussed in Sec. IIHI The 
solutions to the Langevin equation under a dipole approximation can be found through the 
Laplace transform for the charge's motion either parallel or perpendicular to the bound- 
ary plane in Sec. IIVI Thus, velocity fluctuations are obtained in Sec. |VJ The results are 
summarized and discussed in Sec. EO 
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The Lorentz-Heaviside units with h = c = 1 will be adopted unless otherwise noted. The 
metric is rf v = diag(+l, —1, —1, —1). 



II. INFLUENCE FUNCTIONAL AND LANGEVIN EQUATION 

We consider the dynamics of a nonrelativistic particle of charge e interacting with quan- 
tized electromagnetic fields. In the Coulomb gauge, V ■ A = 0, the Lagrangian is expressed 
as 

1 If 
L[q, A T ] = -mq 2 - V(q)-- / d 3 xd 3 y g(x; q)G(x, y)g(y; q) 

1 (1) 



+ /cfx 



i(5 M A T ) 2 +j-A n 

in terms of the transverse components of the gauge potential At, and the position q of the 
charged particle. The instantaneous Coulomb Green's function G(x, y) satisfies the Gauss's 
law. The charge and current densities take the form, respectively, 

g(x; q(t)) = ec^(x - q(t)) , j(ar; q(t)) = e q(t) <j( 3 >(x - q(t)) . (2) 

Let p(t) be the density matrix of the particle-field system, and then it evolves unitarily 
according to 

p(t f ) = U{t f ,U)p{t i )U- 1 {t f ,t i ) (3) 

with U(tf, ti) the time evolution operator of the total system. The nonequilibrium partition 
function can be defined by taking the trace of the density matrix over the particle and field 
variables, 

Z = Tr{U(t f ,U)p(t l )U~ 1 (t f ,t l )} . (4) 

It is convenient to assume that the state of the particle-field at an initial time ti is factorizable 
as p(U) = p e {ti) <S> PA T (ti). The more sophisticated scheme of the density matrix involving 



initial correlations can be found in Ref. 13j. We also assume that the particle initially is in 
a localized state, and thus its density matrix can be expanded by the position eigenstate of 
the eigenvalue q^, 

Pe{U) = \<li,U) (q», *i| • (5) 
The electromagnetic field at the time ti is assumed in thermal equilibrium at temperature 
T = 1/(3, and thus its density operator takes the form 

pA T (t,) = e-^ A T/Tr{ e -^ A T} , (6) 



where Ha t is the Hamiltonian of free vector potentials in the Coulomb gauge. Later, we will 
focus on the case that the initial state of fields is the vacuum state in the zero-temperature 
limit (3 — > oo. The nonequilibrium partition functional can be computed with the help of 
the path integral along the contour in the complex time plane by taking the limits, £j —>■ — oo 
and tf — > +00, and it is given by [ill ] . 



Z 



Q/ ret; 

d 6 f\f I T>q + T>q~ exp 



h 



dt (L e [q+]-L e [q-]) 



where the Lagrangian L e [q] is 

1 



L e [q] = ^mq 2 - V(q) - ^ / (fx d A y g(x; q) G(x, y) g(y; q) . 



(7) 



(8) 



The influence functional T [j + , j ], after tracing out field variables, can be written in terms 
of real-time Green's functions of vector potentials, 

f [j + >r] = ex p ' 



2h 2 



r/W dV [#■(*; q+(t)) <A+ i (x)^(x')>i+(x';q+(t')) 

- j + (x; q + (t) ) <A+* (x) V W) ) 3 J q" (0 ) 

- JT q" (0 ) (A? (s) A* (x') > j+ (x'; q+ (*') ) 

+ j 4 -(x; q -(t)) (^(x) vV))^'; <r(0) 



(9) 



and contains full information about the influence of quantized electromagnetic fields. Here 
the explicit h dependence is restored in the expressions. The Green's functions are defined 
by 



[A+\x)A^{x')) 
{A- T \x)A-\x')) 
[Ap(x)A-\x')) 



[A^A^x')) 6(t - t') + (A^x^A^x)) 0(t' - t) 
;^ T (x')^(x)) 9{t - t') + {A^A^x')) 6{t' - t) 
;4(x')^(x)) =Tr{p AT 4(x')4(x)} , 
;4,(x)4(x')) =Tr{p AT AUx)A{(x')} . 



It is found more convenient to change the variables q + and q to the average and relative 
coordinates, 

q = ^ (q + + q~) , r = q + - q (10) 

The nonequilibrium partition function in terms of the coarse-grained action becomes 

dqj f VqVr exp Us CG [q,r]| , (11) 



where the coarse-grained action Scg reads 

Scalar] = J™ dtr\t)^-mq\t) - W(q(t)) + e 2 V l G[q(t), q(t)] 

-e 2 (^|-^(0V^ l°°dt' G l i[a(t), q (t');t-f} qi(t') 

+ * y /" ^ ^"(0 - ^ v ') - ? m (0v*") g*? [q(t), q(0; * - *1 } + o( 

after being expanded with respected to r. The prime over V- 7 denotes the partial differen- 
tiation over q j (t'). The retarded Green's function and the Hadamard function are defined 
respectively by 

HG%(x-x') = ^(i-t')(^(x),4(a:')]}, (12) 
hG%(x-x') = i({^ T (x),4(x')}>. (13) 

Next, we introduce the auxiliary noise fields £*(£) with the Gaussian distribution, 

V[C(t)] = exp {-? |°° dt |°° dt' [C(t) G%~ 1 [q(t), q(0; * " A ?(0] } , (14) 

and thus the imaginary part of the coarse-grained action can be expressed as a functional 
integration over weighted by the distribution function V[^i{t)]. As a result, we end up 
with 

exp j^Sccfq, r] 

= JvtiPW)]exp{j: Re{S CG hr]}-heJ™ dtr*^j t -q>(t)V^e }. 

The expressions in the squared brackets on the right hand side is defined as the stochastic 
effective action, which consists of the real part of the coarse-grained effective action as well 
as the coupling term of the relative coordinate r l with the stochastic noise 

The Langevin equation is obtained by extremizing the stochastic effective action and 
then setting r l to zero. By doing so, we have ignored intrinsic quantum fluctuations of the 
particle, and that holds as long as the resolution of the measurement on length scales is 
greater than its position uncertainty. The Langevin equation is then given by 

mq 1 + VV(q(t))+eVG[q(t),q(t)]+e 2 ^|-g i (t)V^ 

x dt' Gg[q(f),q(0;*-*'] q J (t') = -he - q l (t)V^j ?(t) (15) 



with the noise-noise correlation functions, 

<f (*)> = 0, (C(t)e (*)) = iGg[q(t),q(0;«-*1 • (16) 

This Langevin equation encompasses fluctuation and dissipation effects on the charge's mo- 
tion from quantized electromagnetic fields via the kernels G 1 ^ and G\ respectively, both 
of which are in turn linked by the fluctuation-dissipation relation [141 ] . The fluctuation- 
dissipation relation is known to play a pivotal role in balancing these two effects in order 
to dynamically stabilize the nonequilibrium evolution of the particle under a fluctuating 
environment. Mathematically, it relates the Fourier transform of the fluctuation kernel G 1 ^ 
to the imaginary part of the retarded kernel G 1 ^ as follows 

Gg[q(t), q(f ); w] = Im {Gg[q(f), q(f )l w] } coth 
In the zero temperature limit, the relation reduces to 

Gg[q(t),q(0;w] = Im{Gg[q(t),q(0;w]} [*M - B{-u)\ . (18) 

It is found that the backreaction kernel functions of electromagnetic fields in the Langevin 
equation ( IT51) appear purely classical due to the fact that the coupling between the charge and 
electromagnetic potentials is linear. The noise-noise correlation functions can in principle 
be computed by taking an appropriate statistical average with the distribution functional 
V[C, l {t)]- It is also seen from Eq. (j!5p that the influence of electromagnetic fields takes the 
form of an integral of the dissipation kernel over the past history of the charge's trajectory, 
as well as a stochastic noise £, which drives the charge into a fluctuating motion. As it 
stands, this is a nonlinear Langevin equation with non-Markovian backreaction, and the 
noise depends in a complicated way on the charge's trajectory because the noise correlation 
function itself is a functional of the trajectory. 

The general solution q % of the Langevin equation can be expressed as its mean value q l h 
and a small deviation 5q l from the mean. Expanding the equation with respect to 5q l and 
then keeping its linear terms, we may decompose the stochastic equation into the equations 
of motion for q l h and Sq l , respectively. The mean trajectory q % h satisfies the homogeneous 
part of the Langevin equation, which describes the purely classical effects. On the other 
hand, the equation for the position fluctuations Sq 1 involves the stochastic noise £\ The 
noise-driven position fluctuations 5q % thus are entirely of quantum origin as seen from an 
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/3u 



(17) 



explicit h dependence in the noise term. The backreaction dissipation effect on the evolution 
of Sq 1 is expected to balance with the effect from the accompanying stochastic noise via a 

ft 

fluctuation-dissipation relation where both effects are of quantum nature [6( . This issue will 
be further studied below. 



III. STOCHASTIC LORENTZ FORCES WITH THE BOUNDARY 



The integro-differential equation (1151) can be cast into a form similar to the Lorentz 
equation. We consider a charged particle moving in the vicinity of a perfectly conducting 
plate. Let the plate be located at the z = plane. Then, the tangential components of 
the electric field E and the normal component of the magnetic field B on the plate surface 
should vanish such that the boundary conditions of the vector potential A are given by 

A = , and A x = A y = , (19) 

leading to 

-af = ° < 20 > 

as the result of the Coulomb gauge. The transverse vector potential At in the z > region 
is given by, 

. , . fd 2 ku f°° dk z 2 f 

+ ct2(k) i k|| (~j s i n ~~ z (— I C()S : 



e *ii-x||-*-* + H.C.,(21) 



where the circumflex identifies unit vectors. The position vector x is the shorthanded no- 
tation of x = (xikz), where xn is the components parallel to the plate. Similarly, the wave 
vector is expressed by k = (k\\,k z ) with uo 2 = + k 2 z . The commutation relations of the 
creation and annihilation operators are satisfied by 

[a A (k),4(k')] = 5 xx ,5(k ll -k\ l )5(k z -k' z ), with A, A' = 1, 2 , (22) 

and are zero otherwise. Then, the retarded Green's function and Hadamard function can be 
explicitly expressed as the sum of the free-space and the boundary- induced contributions, 

Gg(q(*), q(f ); t - t') = i 6(t - t>) J ^ ^ { - e^'^ (23) 
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gio _ _ 2 22 + e ik -[ q W- g (*')]| ( e -M*-*') - C.C. 



lJ 2 



«S(q( t ), q(Oi« " = \ j S I «> t h [f ] { - e-«-*), (24) 

_ _ W _ ^ + 2 ^^.A e ik.[q(0-q(f )] | ( e -Mt-f ) + C.C.) , 

where the position q denotes the location of the mirror image of the original charge at q 
with respect to the boundary at the z = plane, and they are related by q 1 = (5 lJ — 222) q* . 

The integral expression in Eq. (II 5p can be realized in terms of the Lienard-Wiechert 
potential in the Coulomb gauge due to a moving charge, 

^Lw(q) = / dVGg[q(t),s'] (*';q(t')). (25) 



This potential Alw(cl) clearly depends on the past history of the charge's motion 15]. 
Together with the instantaneous Coulomb potential term, the backreaction can be expressed 
in a gauge invariant way to necessarily maintain underlying gauge symmetry, respected by 
the Lagrange we begin with. With the definition of the current density in Eq. (j2J), the 
straightforward algebraic manipulation further shows that the Langevin equation can be 
re-expressed as retarded Lorentz forces and the stochastic components, 



™? ! + V q V(q) = e [J5*(q) + e ijk <f (t)5 fe (q)] - he - q>(t)v£) ?(t) . 



(26) 



The electromagnetic fields are defined by 



E = - J^A LW - V q G , B = V q x A LW . (27) 

Next we write electromagnetic fields in reference to the retarded spacetime coordinates and 
explicitly have them divided into the free-space contribution and the correction out of the 
boundary. Let the free-space part G 1 ^'^ and the boundary correction G^'^ of the retarded 
Green's function be respectively given by 

Gf 0) (q, q'; r) = i 9(t - t>) J |* _L _ e ^-^ + Q Q 

Gg ;(h) (q, q'; r) = -iO(t -if) J ^ 1- - ^ - 222 + 2^^j e ik <^~^ + C.C. 

where q' = q(t') and r = t — t' . The similar decomposition is applied to the instantaneous 
Coulomb potential. Therefore, we end up with 

E = E(°)(q,R) + E( b )(q,R), B= [n x E<°>(q, R) + n x E«(q, R)]^ , (28) 



where 



E (0) (q,R) = e 
E (b) (q,R) = - e 



n — q 



7 2 (1 - q - n) 3 i? 2 
IT q 
7 2 (l-^-n) 3 ^ 2 



+ e 



ret 



n x [(n — q) x q] 
(1 - q- n) 3 R 

n x ](n — q) x q 
(1 - ^-n) 3 E 



ret 



(29) 
(30) 



in which R = q(i) — q(tjt), n = R/i?, and 7 = (1 — q 2 ) -1 / 2 . The retarded time is defined 
as tn = t — R. The variables with the overbar are given by replacing the source point q(tf?) 
with its image position q(t/j) in their definition. The minus sign in the definition of the 
retarded electric field (130]) implies that this field can be interpreted as radiation out of the 
image charge due to the boundary. Thus Eqs. (!28j) - (l30l) are consistent with what would be 
obtained for electromagnetic fields due to a moving charge and its image 151 ]. 

The free-space part of the Lienard-Wiechert potentials involves its associated retarded 
Green's function, which is nonvanishing only for lightlike spacetime intervals. Since the 
worldline of a massive particle is timelike, in the absence of the boundary the charge at 
the present time can not be affected by backreaction from previously emitted radiation due 
to the nonuniform motion of the charge itself Jjj]; however, radiation can backreact on 
the charge at the time when it is just emitted. Thus, the nonlocal form of the free-space 
contribution of the Lienard-Wiechert potential reduces to a purely local effect. Besides it 
may lead to short-distance divergence in the coincidence limit. This ultraviolet divergence 
arises from the assumption that the point-like particle interacts with fields, and must be 
regularized to have a finite and unambiguous result. Then the divergent part is absorbed by 
mass renormalization, while the finite backreaction takes the form of electromagnetic self- 
forces, given by a third-order time derivative of the position. It may bring about issues such 
as the runaway solutions and acausality on the dynamical evolution of a point charge 161 ] . 

Needless to say, the emitted radiation may backscatter off the boundary, and in turn 
affects the charge's motion at a later time. Thus, this backreaction owing to the presence of 
the boundary depends on the past history of the charge's trajectory, leading to the memory 
effect. This non-Markovian processes can also be understood as if radiation was emitted at 
a retarded time from the image charge at q, which is the location of the mirror reflection 
of the point charge at q due to the presence of the conducting plate. The lag corresponds 
to the time delay for radiation to travel from the image counterpart at an earlier time to 
the charge itself at the present time, roughly being equal to the round-trip traveling time of 
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radiation between the plane boundary and the charge. 

In addition, from the interpretation of the stochastic noise and the noise-noise correla- 
tion, we may formally identify —Hg 1 as the stochastic vector potentials A\ such that their 
correlation functions are 

(4(t)>=0, (Al(t)Ai(t')) = hG%[ q (t), q (t'y,t-t']. (31) 

The noise term on the right hand side of Eq ( l26l) can then be thought of as the stochastic 
Lorentz force F s , 

F* = e[El + e l3k qi(t)B k s ] , (32) 

where 

£: = -JU;, and b: = -v;xa;, (33) 

The origin of the stochastic Lorentz force by construction comes from electromagnetic vac- 
uum fluctuations. Here we note that the stochastic electromagnetic fields (E l s ,Bl) involve 
only the transverse components of gauge potentials. In the Coulomb gauge, the instanta- 
neous Coulomb potential, which is determined by the Gauss law, is not a dynamical variable 
when the intrinsic quantum fluctuations of the charge are ignored, and hence it has no corre- 
sponding stochastic component. In the end, the Langevin equation fl26l) for the point charge 
interacting with quantized electromagnetic fields in the presence of the conducting plate can 
be nicely cast into a gauge-invariant form, 

mq l + V^V(q) = e [^(q) + e ijk q 3 (t)B k (q)} + e [&„ + e ijk q> (t)B k ] . (34) 

Then the expressions in the first pair of squared brackets on the right hand side denote 
the retarded electromagnetic forces, while the terms in the second pair are its stochastic 
components, manifested from the quantum fluctuations of electromagnetic fields. 

IV. LANGEVIN EQUATION UNDER DIPOLE APPROXIMATION AND ITS SO- 
LUTION 

The nonlinear, non-Markovian Langevin equation is far too complicated to study fur- 
ther without any approximation. We will consider that a charged particle undergoes the 
harmonic motion, and assume that the amplitude of oscillation is sufficiently small. The 
appropriate approximation for a non-relativistic motion will be the dipole approximation. 
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This approximation amounts to considering the backreaction solely from electric fields, and 
linearizing the Langevin equation in such a way that the equation of motion remains non- 
Markovian. Since the presence of the boundary will anisotropically modify the effects of 
electromagnetic fields on the evolution of the charge, we will consider the motion in two 
different directions, that is, either parallel or perpendicular to the plane boundary. 

A. parallel motion 

When the charged particle moves parallel to the boundary, say, in the x direction, let the 
equilibrium point be located at the coordinates (x, y, z) = (0, 0, Zq). The linearized Langevin 
equation for the motion displaced from its equilibrium position reduces to 

mf(t) + d 2 x V(z ) q x (t) + e 2 f di! g x R [z , z 0] t - t'\ q x (t') = -Ke?(t) , (35) 

J o 

from Eqs. ( |T5|) . ( |23l) . and ( 12^1) under the dipole approximation. Here q % denotes the dis- 
placement from the equilibrium point. The xx component of the retarded Green's function 
in the dipole approximation is denoted by g R [z ,z ;t — £'], and can be written in terms of 
the spectral density p x as 

r°° dk 

g x R [z , z ; t-t'] = -6{t - t') — p x {z , z ; k) sva[k{t - t')} , (36) 

Jo K 

where the spectral density is given by 

1 sin(2A;2 ) cos(2A;z ) sm{2kzo) , . 

3 ~ 2(2kz ) ~ 2(2kz ) 2 + 2{2kz f \ ' ^ ' 

The instantaneous Coulomb potential in this case is found to have no effect on the motion in 
the x direction, but will establish a static attraction force between the charge and its image in 
the z direction. Hence, the applied potential V is assumed to have an additional component, 
other than the harmonic potential, to counteract the static force so that the motion of the 
charge remains on the z = z$ plane. In addition, the motion has been assumed to start at 
t = 0, and its initial conditions are chosen to be q x (t < 0) = q x (0) and q x (t < 0) = 0. In 
general, it is insufficient for a non-Markovian equation to have a unique solution if a finite 
number of initial conditions are specified. Thus in this case, the initial conditions are given 
over half of the real axis t < 0, although later it will become clear that the requirement 
on specifying initial conditions can be less stringent in the current case. Physically, these 
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p x (z ,z ; k) = 

7T 



initial conditions can be achieved by applying an appropriate external constraint to hold 
the particle at rest at the position q x (0) for t < 0. Then the applied potential is suddenly 
switched off to V(q) at the time t = so that the dynamics of the particle at later times 
follows the Langevin equation. The accompanying noise-noise correlation functions due to 
quantum fluctuations of E x fields are given by 

(C(t)) = 0, (f (W)> = IgU^ZoM-t'], (38) 

where 

f°° dk 

g x H [z , z ;t-t'] = - — p x {z , z Q ; k) cos[k{t - t')} . (39) 
Jo 271 

The noise kernel g x H is related to the dissipation kernel g x R via a fluctuation-dissipation 
relation under the dipole approximation. 

Carrying out the integration over k makes the Langevin equation physically more trans- 
parent, 

q x (t) = -hei x (t) + pq x (t) (40) 

q x (t-T) q x (t-T) q x (t-T)~ 
(2* ) + (2^ ) 2 + (2z ) 3 . ' 

with the time delay T = 2zq/c and the renormalized mass m r = m + e 2 A/37r 2 . Here for t > 
the charged particle is set into harmonic oscillation with frequency Uq by specifying the 
proper potential V(q). The ultraviolet divergence arises due to summing up backreaction 
effects from the free-space contribution over all energy scales of fields in the coincidence 
limit. This type of divergence is often seen in quantum theory as well as classical theory 
when an infinite number of degrees of freedom are involved. The energy cutoff A is then 
introduced to regularize the integral. The cutoff scale can in principle be chosen to be the 
inverse of the particle's classical radius. However, this energy scale is already well beyond 
the regime of the validity of a nonrelativistic, single-particle description for the Langevin 
equation. Thus, instead, the more sensible energy cutoff scale will be the inverse of the 
width of the charge's wavefunction. It essentially quantifies the intrinsic uncertainty on the 
charged particle. The divergence is then absorbed into mass renormalization as shown above. 
The finite backreaction effect is given by the third-order time derivative of the position. 

The boundary-induced backreaction in the squared brackets on the right hand side of 
Eq. fj40|) reveals the non-Markovian nature in the argument of each individual term. The 
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m r q x (t) + 



m r (jj Q + 



4tt (22 )3 




assumption of small amplitude for the charged oscillator simplifies a general non-Markovian 
integro-differential equation into an ordinary differential equation with a fixed time delay. 
This is expected due to the fact that by small amplitude, we mean the variation of the 
position for the charge oscillator is much smaller than its distance to the plate Zq. It implies 
that the time variation owing to the oscillation amplitude is ignorable compared to the time 
2zq/c needed for the radiation to have a round trip between the charge and the boundary. 
Therefore, the dominant contribution of the time delay will be just given by a fixed constant 
T = 2z /c. The term with q x (t — T) is derived from the acceleration field of electric fields 
while the rest two terms come from the velocity field. The time difference t — T indicates 
that the backreaction effect due to the boundary on the charge at the time t depends on the 
charge's dynamics at an earlier time t — T. The time delay 2z /c may also be thought of as 
the traveling time for radiation emitted from the image charge at an earlier time to reach the 
charge at a later time. The memory effect is thus described by such a time delay differential 
equation. This equation can also be obtained from Eq. ( |26l) by expanding the expression of 
the retarded Lorentz force to the linear terms in q x , q x and q x in the non-relativistic limit. 

In order to find the solutions, it apparently requires the knowledge of q x at earlier times. 
Since the Langevin equation under consideration reduces to an ordinary differential equation 
with a fixed time delay T, the specification of the initial conditions can be relaxed to the 
conditions 



in the interval — T < t < 0, instead of the whole negative real axis of t. Then the equation 
can be solved by means of the Laplace transform. From Eq. the Laplace transformed 
Langevin equation reads 



m r [s 2 r(s) - s q x (0) + + e 2 ^[z , z ; s] [s 2 F(s) - s q x (0)] = -fieff (a) , (42) 



where the initial conditions (T4T]) have been applied. The Laplace transformation of a function 
fit) is defined by 



<f(-T <t<0) =<f(0), 



q x (-T <t<0)=0 



(41) 




t > 



while the inverse Laplace transform of f(s) is 
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where the Bromwich contour C counterclockwisely encompass all singularities of f(s) on the 
complex s plane. Since ultraviolet divergence has been absorbed into mass renormalization, 
in order to avoid being cluttered with the notations all variables will be assumed renormalized 
from now on unless specified otherwise. The renormalized kernel g%[zo, z$\ s] is given by 

dh T , T * k A 
— P [zo,z Q ; k) — 



The solution to the Laplace transformed equation is readily obtained as 

s{l+ e i9 x R [^z ;s]}q%0)-^(s) 
§*(s) = -± ' 



(43) 



where the T, x (s) kernel is 



S 2 + UJ 2 + Z x (s) 



E x (s) = — s g%[zo, z ; s 
m 



(44) 



(45) 



Thus, the solution in the time domain can be given by the inverse Laplace transformation 
of q x (s), and it reads in terms of the mean trajectory and its deviation from the mean, 



q x (t) 



,2 rt 



K x (t) + - / dt' g x R [z , zo]t- t'] K x (t') 



eh r l 

<f(0)-— / dt' K x (t-t')C(t'), (46) 



m 



where the kernel K x (t) is the inverse Laplace transform of K x (s) = [s 2 + uJq + S x (s)] 1 

1 



K x (t) 



ds 



(47) 



lc s 2 + uj 2 + E x (s) 

with K x (0) = and K x (0) = 1. It is seen from Eqs. (@3} and (SSD that the kernel E x {s) 
has a branch-cut along the imaginary s axis, so the discontinuity of the kernel E^s) over 
the branch-cut can be made explicitly by letting s = iu + (J 1 * 1 . In terms of uj, the kernel Y* x 
becomes 

S x (s = i u + C^) = Re T, x (uj) ± i Im T, x (uj) , 

where 



(48) 



Re Tj x (uj) = — uj 1 
m 



UJ' 



A 

3^2 



ImSj-fw) = sgn(u;)u; 2 / dk k p x (zo, Zo; k) 8(k 2 — uj 2 ) 

m Jo 



Carrying out the integrals yields 



Re £ T (u;) 



■ UJ 



I cos(2u;z ) sm(2uz ) cos(2ujz ) 



_(2ujz ) 3 (2uz ) (2ujz ) 2 (2ujz ) 3 _ 



(49) 
(50) 

(51) 
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FIG. 1: (a) On the first Riemann sheet, the branch-cut lies along the real axis of u and a runaway 
pole sits at the negative imaginary axis. All the other poles stay on the second Riemann sheet. 
Thus, the Bromwich contour reduces to two straight lines parallel but infinitesimally close to the 
cut from both sides if the contribution of the runaway pole is discarded, (b) The resonance modes 
are found asymptotically in the complex uj plane along the curve C p . 



JmEJu) 



A-nm 



sgn(u;j uj" 



2 sin(2co>,2o) cos(2c<jzo) sin(2u;zo) 

3 (2uzq) (2ujz ) 2 + (2tuz ) 3 



(52) 



In general, the branch-cut of the S a (s) kernel lies within the intervals (u;th,A) and 
(—A, — Uth), where A is the energy cutoff and uth is the threshold energy above which the 
process of particle creation is possible [17|. However, for quantized electromagnetic fields, 
its threshold energy vanishes due to masslessness of the photon, i.e. u> t h = 0. Thus the 
branch-cut will extend along the imaginary s axis from — iA to iA. 

The locations of the poles on l/(s 2 + uJq + H x (s)) can be found by solving the following 
equation, 

s 2 + ul + E x (s) = 0. (53) 

This pole equation turns out to be a transcendental function of s due to the presence of the 
boundary. It may have an infinite number of solutions. This in turn implies that the integro- 
differential equation can be viewed as an infinite order of the ordinary differential equation. 
Thus an infinite number of initial conditions are needed 18]. For the weak coupling constant 
e 2 , perturbatively, we may let s = i{ujq + 5uS) + Sj, where 5uj — i5^j is the shift of the pole 
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owing to interaction with electromagnetic fields. It is then plugged into the equation ([53]) 
to find Sou and £7. We assume £7 > for the moment. To the order e 2 , the perturbative 
solution is found to be 8uj = ReT, x (u ) / (2uo) , and £7 = -ImE3.(w )/(2w ). However, since 
from Eq. (151Z1) . ImSj.(w) > for all positive values of lu, it leads to £7 < 0, inconsistent with 
the assumption we made earlier. A similar contradictory result will be found through the 
same argument by assuming £7 < 0. As a result, the poles do not exist in the first Riemann 
sheet, and may appear in the second sheet. This is the known fact for the resonance {20]. In 
fact, there are still an infinite number of poles lying on the upper half of the second Riemann 
sheet. However, one exception is the runaway pole, which resides on the positive real axis 
of the complex s plane. Since the kernel S x (s) is real when evaluated at the runaway pole, 
this pole stays on the first Riemann sheet. It corresponds to either a runaway motion of the 
charge or the acausal evolution due to preacceleration. Since it is viewed as unphysical if the 
asymptotically bounded motion is considered we then discard the contributions of the 
runaway-type poles to the inverse Laplace transformation ( 1471) . Accordingly, the contour 
on the first Rienmann sheet can be deformed to be parallel and infinitesimally close to the 
branch-cut as shown in Fig. [U The real-time solution of K x (t) can be written as 

KM) = / 9 sin (ojt) , t > . 54 

J„ th=0 2tt [0,2 _ w g _ Re S ^)] 2 + [ Im s^)] 2 

Similarly, the time derivative of the K x kernel is given by 
• . . [°° duo AojImEJoj) 

The integrand in a weak coupling limit shows a Breit-Wigner feature. The width of the 
resonance is related to the imaginary part of the kernel H x , and its peak is located at the 
resonance frequency Q x to be determined later. For the sufficiently late times, t ^> T, the 
time derivative of K x (t) can be approximated by taking into account only the contribution 
of the resonance mode, and is given by 

K x (t) ~ Z x cos(n x t + a x ) e" r ** , t > , (56) 

where the resonance frequency Q x and the decay constant T x are given by 

ReS g (^ ) ImZ x (n x ) 
il x ~u -\ — , L X ~Z X — , (57) 



17 



respectively. In addition, the phase shift a x and Z x are 



9ReS x (O x )l 1 3imE x (Q 



^-Z x —— . (58) 



The long-time dynamics of the kernel K{t) can also be investigated by examining the low- 
frequency behavior of the integrand in Eq. (1551) along the branch-cut in the neighborhood 
of threshold energy. When the resonance peak is far away from threshold energy, cu t h ^ Q, 
the imaginary part of the kernel E generally behaves like Im E oc (u — w t h) n as u approaches 
Wth- For nonzero threshold energy, uj^ 7^ 0, the integration over u in Eq. (1551) leads to 

f°° 1 n+1 
K(t) oc / dcj (u - u)th) n cos(ut) = — -j- cos(u;tht H ^— 7r) , (59) 

which results in a power-law decay for its late-time dynamics [19] . However, for electromag- 
netic fields with vanishing threshold energy, infrared photons can be generated even with 
an infinitesimally small amount of energy. The backreaction is expected to damp out the 
motion of the particle more effectively than that of the massive field. Thus as shown in 
Eq. (152p . because ImE is proportional to u 3 for small u, we find 

POO 

K(t) oc / duj uj A cos(ujt) = , (60) 
Jo 

after the proper regularization. This result can also be obtained from Eq. ( 159|) . The con- 
clusion holds true for all higher-order time derivatives of K(t). The kernel K(t) thus decay 
faster than the power law, and its time derivative K(t) at asymptotical times is described 
by an exponential decay given by Eq. (1561) . 

To account for the dynamics of K(t) at all times, we need the knowledge of all solutions 
to the relevant transcendental equation u 2 — uJq — ReE^u;) — ilmE^u;) = 0. It is highly 
improbable to analytically locate all of them; nonetheless, when \u\is much larger than unity, 
with the help of Eqs. (15T|) and (|52|) . this equation may be transformed into a exponential 
polynomial such that the asymptotic solutions can be found. We decompose the solution u 
into its real part u and imaginary part v, i.e. u = u + iv. The coupled equations of u and v 
can be given by 



In Vu 2 + v 2 — 2z v = In - — , 

Az 

tan" 1 - + 2z u = (2n + l/2)vr, neZ. 
u 
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Thus asymptotically, the solutions are evenly distributed in the u sense along the curve 
In \/u 2 + v 2 — 2zqv = ln(3/4zo) on the u — v plane as shown in Fig. [TJ With \u\ \v\, both 
u and v are given by 

u = ^(2n + l/2)7r , H>1, 



v = — — In h In m > . 

It can be seen that v is always positive so that all the asymptotical solutions are on the 
upper half side of the complex u plane. These modes will contribute to the motion via the 
time evolution factor 

exp(iuot) oc n~T exp(i 2mr — ) , \n\ ^> 1 . (61) 

Their contributions are thus negligible for the dynamics of K x (t) at times t > T. These 
relatively high frequency modes are significant only in the very early stage of the evolution. 
Therefore, the resonance mode, which peaks at Q with width T, have the dominant effect for 
the late-time dynamics of the non-Markovian Langevin equation in a weak coupling limit. 
The memory effect or the effect of the time-delay will register in the parameters such as the 
resonance frequency Q, the decay constant T, and so on. 

The Laplace transform method seems too complicated to properly take all asymptotical 
resonance modes into consideration in order to deal with the very early-time behavior of 
the charged oscillator. It would be more straightforward to apply the iteration techniques. 
The idea goes as follows. From the initial conditions in the interval — T < t < 0, we may 
find the solution to Eq. ( 1401) for the times < t < T. Because in this time interval, the 
retardation effect does not settle in yet, the equation of motion essentially describes the 
Markovian motion due to the damped oscillator with a driving source. Next, this solution is 
then plugged into the retardation terms on the right hand side of Eq. fT40l) to iteratively find 
the solution for the time interval T < t < 2T. We may continue with these procedures, but 
the secular terms are expected to appear after several iterations. The method of dynamical 
renormalization group can then be implemented by resumming these secular terms to address 
more on the non-Markovian nature of the full-time dynamics of the solution 19]. 



19 



B. perpendicular motion 



We now consider a charged particle moving in the z direction perpendicular to the plane 
boundary. Contrary to the previous case, the modified instantaneous Coulomb term due to 
the boundary will make a shift to the oscillation frequency of the motion. Besides, it causes 
a static attraction between the charge and its image, so an applied potential is introduced to 
counteract the attractive Coulomb force. The corresponding linearized Langevin equation 
for the displacement from the equilibrium position is obtained from Eqs. (iTBI) . (123"]) . and ([21 
as 



mq z (t) + 



d 2 z V(z ) 



4 



4tt (2z ) 3 \ 



+ 7-7^ + e ' / dt ' 9rW z ; t - t'} q z {t' 
47T {2z ) z Jo 

= -het(t), 



(62) 



where the zz component retarded Green's function g z R can be expressed in terms of the 
spectral density, 

f°° dk 

g z R [z , z 0] t-t'} = -6(t -t') — p z {z , z 0] k) sm[k{t - t')\ , 

Jo 7r 



with 



p z (z ,z ]u) 



k 

71 



1 cos(2fcz ) sin(2kz 
3~ {2kz ) 2 + {2kz ) 3 _ 



The noise-noise correlation functions are then obtained as 



<**(«)> =0, 
9m z 0i z o'i t — t] — — 



1 



dk 



(Z z (tmt')) = -gUz ,z ;t-t'} 
p z (zo,Zo;k) cos[k(t-t')] . 



(63) 

(64) 
(65) 



, 2vr 

As discussed before, a fluctuation-dissipation relation is obeyed. 

After mass renormalization, a similar time-delay differential equation for the charged 
oscillator moving in the z direction in the presence of the boundary is given by 



m r q z (t) + 



m r u 



4tt (2z )3J 



tf(t) = -hei*{t) + — q'(t) 
bn 



+ 



2e^ 

47T 



q z {t-T) q z {t-T) 



(2^ 



+ 



(2z f 



The applied potential V(q) has been turned on to drive the charged particle into a harmonic 
motion of frequency ujq in the z direction. In addition, the retarded Green's function also 
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has the contribution to the frequency shift. This equation explicitly shows the memory 
effect due to the boundary. In particular, all backreaction terms arising from the boundary 
come from the velocity field. The contribution of the acceleration field vanishes under the 
dipole approximation since the relative direction from the image charge to the charge itself is 
perpendicular to its motion. Then, the real-time solution q z to the inhomogeneous stochastic 
equation (|62|) takes a similar form to Eq. (H6|) . The time evolution of the kernel K z (t) from 
the inverse Laplace transform is obtained as 



K z {t) ~ Z z cos(f2 z t + a z ) e 



-r z i 



t > 



(67) 



with the resonance frequency Q z and the decay constant T z given by 



n z ~ ul 



+ 



ReE z (^ ) 



Im £ z (f2 2 



(68) 



4irm r (2z ) 3 ' 2cv ~* * 2tt z 

respectively where the frequency shift due to the boundary corrections of the instantaneous 
Coulomb potential has been taken into account. Moreover, the phase shift a z and Z z are 
found to be 



1 - 



dReZ z {tt z 



n -1 



a z ~ Z 



«9Im£ z (a 



(69) 



Both the mass and the real part of the £ z (a>) kernel have been renormalized. The corre- 
sponding £ z (s) kernel is 



m 



s 9r[ z o-i z o', s ] 



e 
m 



dk s A 

— fcp [Zq, z ; k) — — — — - 

7r s z — k z 3n z 



(70) 



and the real and imaginary parts of S z (s) in the vicinity of the branch-cut are given, re- 
spectively, by 



Re TiJu) 



Im Ti z (cu) 



e 2 



2 sin(2cjzo) 2 cos(2a;zo) 



(2uz ) 3 (2uz y 



(2a; zq) 



sgn(u;J ijj 



2 2 cos(2u;zo) , 2 sin(2a;2;o) 

3 " 



(71) 
(72) 



(2^ ) 2 (2cuz ) 3 

The presence of the boundary apparently modifies the behavior of the charged oscillator 
in an anisotropic way. This characteristic is especially noticeable near the boundary where 
the electric fields parallel to the plate vanish, but their normal components become doubled, 
compared with their counterparts without the boundary. This anisotropic feature is encoded 
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in the spectral density. Thus, the quantities that can be expressed with the spectral density 
should share the same feature. As a result, the decay constant T for motion parallel to the 
boundary turns out to be smaller than in the perpendicular case with a similar configuration. 
Additionally, the effect of the stochastic noise on the oscillator ought to be much weaker in 
the parallel case than the perpendicular one. 

Next section will be devoted to studying the time evolution of velocity fluctuations of the 
charged oscillator due to both anisotropically modified vacuum fluctuations by the boundary 
and associated dissipative backreaction from electromagnetic fields. 



V. VELOCITY FLUCTUATIONS 



When the charged particle couples to quantized electromagnetic fields, the stochastic 
Lorentz force, manifested from vacuum fluctuation of fields, drives the particle into a fluctu- 
ating trajectory in analogy to the Brownian motion. Thus, it is of interest to study velocity 
fluctuations of the charged oscillator to see how they are affected by the boundary and are 
asymptotically saturated as a result of the fluctuation-dissipation relation. We will compute 
them in the interval l/fi<Ct<l/r for the linearly growing regime in which backreaction 
dissipation can be ignored, as well as in the interval of 1/T <C t for the saturation regime 
where the effects of fluctuations and dissipation come into balance. 

As the particle starts to move at t — 0, its velocity fluctuations at time t in the direction 
i driven by electromagnetic vacuum fluctuations are given by 

(Aiftf)) = ^ J* dt' ^ dt" [g'nizo, *b5 t' - O] KS") > (73) 

where g l H is the Hadamard function of vector potentials in the dipole approximation. We 
have implicitly assumed that (Aw 2 ) vanishes initially. The charged oscillator is allowed to 
move either parallel or perpendicular to the boundary, and thus the relevant Hadamard 
function takes the form of Eq. ( |39l) or ( l65l) . The function Ki(t) is the kernel of the equation 
of motion, and its time derivative is approximated by Eq. (15 61) or (!6"TI) . This approximate 
solution holds for t > 1 /0» > 2zq by ignoring the contributions from high-frequency modes. 
Then, it may give rise to errorr in computing the integral f!73l) for the time regime t < 2zq. 
In fact, one can argue that the error for a large distance z can be trivially neglected because 
the boundary correction is negligible. Moreover, although the presence of the boundary will 
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result in significant contributions for small Zq, the time interval < t < 2zq over which 
integrations in Eq. f!73|) are performed, are also small. Only a minor error is introduced. 



In terms of the spectral density, the velocity fluctuations becomes 

,2 rt rt roo 



p 2 ft rt roo 11 

(Av 2 (t)) = / dt' dt" KiWKtf') / — k 2 p\z , z ; k) cos[k{t' - t")} 

171 Jo Jo Jo ^ n 

p 2 f°° dk 

= -J ^- k 2 p\z Q , z q ; k) I(t; k) , (74) 

m Jq Z7t 

where 

I{t-k)= [ dt' [ dt" ki{t')ki(t") cos[k{t' - t")] . (75) 
Jo Jo 

In writing so, the integrand of the /c-integral is factorized into a product of the spectral 
density p l (zo, Zq] k) and the function I(t; k), which solely depend on the distance zq and time 
t, respectively. The full expression of the function I is too intricate to be of any help, but 
in the time intervals 1/Qi <C t <C 1/Tj and t ^> 1/Tj, it can be reduced greatly, 

1. <C t <C 1/Tf linearly growing regime, with T^t set to zero, 

I (t,k) = I(—<t<^;k) 

2(k 2 + tt 2 ) + (k 2 - tt 2 ) cos(2ai) + (k 2 - tt 2 ) cos(2^ + 2a, t ) 
- (k + Vti) 2 cos[(/c - tti)t\ — (k — tt { ) 2 cos[(k + tt^t] 



2(k 2 -n 2 ) 2 



- (k - ttf) cos(fi;t -kt + 2cti) - (k - ttf) cos(^t + kt + 2a { 

2. Fit 3> 1: saturation regime, with Tit set to infinity, 

2 (k 2 + rf ) cos 2 ttj - 2QiTi cos a» sin a, + tt 2 sin 2 a» 
loo(k) = 1 [t —> oo; k) — Z i 



[(k - a) 2 + r 2 ] [(k + + r 2 ] 

Notice that the function I has a Breit-Wigner feature in the k space with the resonance 
peak at about tti and its width being approximately of order 7r/t at early times or Tj for 
the late-time regime. The narrow width may result in a prominent peak. On the other 
hand, the spectral density p % reveals the oscillatory behavior on the scales 1/zq in the k 
space, which will lead to heavy cancelation as the distance zq is sufficiently large so that 
the boundary effect becomes insignificant. Thus, the result of the integration (1741) relies on 
two competing scales 1/zq, and ir/t or IY Additionally, the behavior of the integrand in 
Eq. ( 1741 increases linear in k when k is sufficiently large. Therefore, velocity fluctuations 
after doing integration over k will have a quadratic divergence, which has to be regularized 
by introducing a cutoff. Zi is approximately equal to unity in the weak coupling limit. 
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A. Growing regime 



Here we study the motion of the charged oscillator at the early stage for <C t C 1 /I\ 
when dissipation backreaction can be ignored. Velocity fluctuations thus mainly result 
from the stochastic noise. For t ^> 2z when the retardation effect out of the boundary 
becomes effective, we may expect that the spectral density varies relatively slowly with k in 
comparison with the function / around the resonance peak. 

Therefore, the integration (174|) can be approximated by pulling the spectral density p % 
out of the integral and evaluating it with k ~ f2j in the neighborhood of the peak. With 
this in mind, we then rewrite Eq. (1T41) in terms of the dimensionless parameters denoted as 
7i = Ti/Qi, z 0i = f^o, y = k/Qi, and Tj = Q{t, and obtain 

e 2 tf r dy 
'o 



<A^ 2 (t)> 



2n 



V p l (z 0i ,z 0i ;y)l(Ti;y) , 



(76) 



where 



f 



p*(zb», zof, y) = ^P*(^o, z ; k) 



l{T i ;y) = Q 2 i Io(t; k). 

Since the function 1 demonstrates a sharp peak about y = 1 with the width of order 7r/rj, 
as long as Zqi/t-i <C 1, the velocity fluctuations can be approximately given by 



■ 2 '^ ~ '^ 2 '^ e2nj p\z 0u zo^) rp-z(n;v) 



(Avt(t)) ~ (Avf(t)) div . - 



2tt 



where 



— T{ji\ y) « — + terms oscillatory with time 
o 2tt 4 



(77) 



(78) 



leads to a contribution linearly growing in time. The high-frequency contributions 
(Avf(t))div. need regularization by inserting a convergent factor e~ yei into the integrand. 
Here 6j is the dimensionless short-distance cutoff. Then, it ends up with 



2 + cos(2aj) + cos(2^t + 2a { , 



1Q-K 2 m 2 



3 + cos(2fi i) 



3A dB 
4 

3A dB 



+ 



7? > 



4 

4 , 1 

x direction 
z direction . 



x direction 
z direction . 



(79) 



Notice that to obtain the last expression, the dimensionless short- distance cutoff e, is chosen 
as a product of the resonance frequency and the de Broglie wavelength of the charge AdB, 
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that is, €i = fijAdB- The cutoff dependence terms are known to arise from the unbounded 
Minkowski vacuum fluctuations, experienced by the charged oscillator. They weakly depend 
on the modified oscillatory motion that results from the coupling to vacuum fluctuations 
in the presence of boundary We will make further approximations by $7j ~ Qq, oti ~ 0, 
and Zi ~ 1, ignoring higher order effect on z . For motion perpendicular to the boundary, 
the dominant zq dependence of high frequency contributions is thus the 1/z 2 term, and 
should be small as compared with the cut-off dependence term in accordance with the 
classical assumption on the particle. However, the corresponding term is missing for motion 
parallel to the boundary. This can be understood as the consequence that there exists no 
corresponding boundary for the motion in the x direction; thus no length scale is introduced 
in this direction. 

Now we consider the ratio of the cut-off dependence term to the term linear in time, ob- 
tained form the above expression. With the help of the explicit expressions of the relaxation 
constants Y h given by Eqs. (I52l). (jBTjl . (j68l). and (j72|) . 

(80) 
(81) 



ryoo) = -— n 



< - „ 2 2 sin(2fio-2 ; o) cos(2f2oZo) t sin(2f2o^o) 
Tcm 



87rm 



3 (2Q z ) {2Q zo) 2 (21Vo) 3 
2 2 cos(2Q z ) ^ 2sin(2fi <2o) 



3 (2fi ^o) 2 (2fi ^o) 3 
it can be shown that the ratio is of the order 

11, 1 Is 



Apparently, the high-frequency contributions can be dominant for the sufficiently low oscil- 
lation frequency, or at very early stage of motion 12]. At later times, e -2 f2o -1 <C t, they 
then become insignificant as it can be further seen from the numerical result to be discussed 
later 

Therefore, the velocity fluctuations are found to grow linearly with time. The growing 
rate is related to the relaxation constant Tj, evaluated at the resonance frequency Qq. Back- 
reaction dissipation will set in at later times, and will asymptotically counteract the effects 
of fluctuations. 
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B. Saturation regime 



Next we will investigate the behavior of velocity fluctuations at much later times t ^> 1 /Tj 
by incorporating backreaction dissipation. In this regime, they are described by 



<A^ 2 > 



e 2 a 



mr 



7T y 2 p l (z 0i ,zoi-y)l(y) , 

ZTT 



(82) 



where 



l(y) = tf /oo(fc). 



In general, since Eq. (1821) has a quadratic divergence, a regulator of the form e~ yti will be 
introduced to regularize the integral. The influence from the boundary can be significant 
when the dimensionless distance of the oscillator from the boundary z j is much smaller than 
unity. 

Thus, we will expand the result of velocity fluctuations in terms of small zoi, after we 
plug in the dependence of the parameters ji, a>i, and p l for the motion in the direction 
i, which are given by, 



P x (z x,z 0x ;y) 
p z (z z,z 0z ;y) 
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= -H 

71 

cos(22o; 



1 _ sm(2z 0x ) 
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3 



+ 



2(2z 0x ) 
cos(2z 0z ) 

sm(2z 0x 
2(2%r) 



+ 



+ 



cos(2z 0x ) t sin (2z 0x 

2(2Zoa;) 2 

sin(2z ^) 
(2z 02 ) 3 



2(2^o x ) 



1 + 



sin (2-2, 



{2z 0z 



and the dimensionless relaxation constant can be obtained from Eqs. (|80|) and (IHTI) as j$ = 
—7if ci p l with f ci = QiT c . The parameter r c = e 2 /(47rm) is the particle's classical radius. 
Here, another small parameter f C j, characterizing the distance r c over the oscillation time 
scales due to a nonrelativistic motion, will also be used to extract the dominant contributions. 

The saturated value of velocity fluctuations can be found analytically. For the motion 
parallel to the boundary, it is given by 
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For the perpendicular case, we find 
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where r y e is Euler's constant, with numerical value ~ 0.577216. The last lines of Eqs fl83|) 
and (184p are further approximated by Zj ~ 1, f2j ~ fin, ignoring higher-order zn-dependent 
terms where z = fin^o, f c = finr c , and e = QoXdB', hence typically z > e ^> f c . The 
parameter fi is the renormalized oscillation frequency with the frequency shift due to the 
interaction with fields. Here for simplicity, the same Q is chosen for the motion in both 
directions because the anisotropy it introduces is of next order in e 2 . In general, the re- 
sults ( J83l) and ( |84l) should show strong dependence on the relative orientation between the 
boundary and the direction of the motion for small zo, and thus are anisotropic. It can be 
understood as a result of vacuum fluctuations of electromagnetic fields under the influence 
of the boundary ll( . The enhancement in velocity fluctuations due to the presence of the 
boundary in the direction normal to the conducting plate arises from large induced E-field 
fluctuations near the boundary, in comparison with the fluctuations parallel to the plate 12]. 
On the other hand, the term depending on Qq results from particle's motion, and is found 
to be boundary-independent. The physics behind this features can be explored in a more 
general context as follows. Starting from Eq. ( |82|) . for small zoi, and f C j, the contribution 
from the resonance peak is given by approximating the integral over y with the spectral 
function evaluated at the dimensionless resonance frequency, y ~ 1, 
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where again the approximation Qi ~ Qq has been applied to obtain the last expression. Ap- 
parently, the result superficially depends on the spectral function p % \ however, the integral to 
the right of the spectral function yields an expression proportional to l/7j. Since the dimen- 
sionless relaxation constant 7, can be related to the spectral density by — nf ci p l (zoi, Zen; 1), 
the boundary dependence due to the spectral density p 1 is canceled. Thus, the result does 
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not depend on the relaxation constant, hence independent of the orientation. Recall that 
the function p 1 is proportional to the rate at which velocity fluctuations grow, while the 
time scales to reach the saturated regime are given by 1/7;. Thus it implies that although 
the velocity fluctuations increases at different rates at early times, they will take different 
amount of time to reach saturation. In the end, the saturated values are independent of the 
direction of the motion, and are only determined by the oscillation frequency of the charged 
particle. This cancelation of the boundary dependence, in fact, can be understood as rem- 
iniscence of the fluctuations-dissipation relation since the function p 1 comes from the noise 
kernel g l H while the dimensionless relaxation constant ji is inherited from the dissipation 
kernel g R . 

As stated above, the integrand of Eq. (|82|) grows linearly for large values of y, so the 
contribution of this portion is given by 



where e = VloXdB, and the lower limit a of the integral is chosen to be greater than unity but 
otherwise arbitrary in order to exclude the integration region around the resonance peak. 
As long as z 0i a <C 1, the integral is pretty much independent of the choice of a. Clearly, 
the integral flHBI) is boundary- dependent and thus anisotropic. The quadratic dependence 
on the cutoff can be characterized by the width of the particle's wave wavefunction. When 
the spatial extension of the wave function is smaller, it tends to probe the finer structure of 
fluctuations, which in turn causes larger variation in particle's velocity. 

Some comments are in order. In the limit Qq — > 0, the motion of the charged particle 
becomes insignificant and the dissipation effect is ignored. Hence, it is of interest to compare 
our results with the earlier study [12]. They are consistent as t — > 00, apart from a factor 
1/2 due to an average over a period in the oscillatory motion. In this limit, the backreaction 
dissipation can be safely ignored because the relaxation constant Tj with f2g dependence, is 
vanishingly small for such a slow motion. Nonetheless, for a finite value of Qq, the nonuniform 
motion results in a dissipative effect, which damps out the particle's motion. 
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C. Numerical results 



The typical behavior of time variation of velocity fluctuations is shown in Fig. [21 where 
the results are generated numerically. The value of the cutoff is set by the de Broglie 
wavelength, X dB ~ h/{mv), so the ratio of the charge's classical radius to its de Broglie 
wavelength is r c /AdB ~ {e 2 /he) {v/c) ~ 10~ 4 with v/c ~ 10~ 2 for nonrelativistic motion. In 
Fig. [21 we choose z = z fl Q /c ~ 1, for example, when the distance to the plate z is of the 
order /im, the oscillation frequency is given by Qq ~ 10 14 s _1 below the plasma frequency 
of the plate 10 16 s _1 [l^. Now for an arbitrary charged particle, without loss of generality, 
we may take appropriate values of the parameters to numerically compute full evolution of 
velocity fluctuations. Here we let m ~ 10 3 m e where m e is the mass of an electron, leading 
to the value of e = AdB^o/c ~ 0.1 and f c = r c Qo/c = 5 x 10 -4 . The horizontal axis is 
the dimensionless time t = Q t, and the vertical axis is velocity fluctuations, normalized 
by the parameter h~Qo/2m. It is seen that at early times, the velocity fluctuations increase 
linearly at different rates for motions parallel and perpendicular to the conducting plate, 
and then saturate to the same value at late times. In this case the cutoff-dependent terms 
is negligible except for very early times t < e~ 2 Q$ 1 ~ lO 2 ^ 1 , which is vanishingly short 
when compared with the relaxation time scales for motion in either direction, about the 
order of 1/T ~ lO 4 ^ 1 . Thus, their contribution can not possibly be seen in Fig. [21 



VI. DISCUSSIONS AND CONCLUDING REMARKS 

The anisotropic behavior of velocity fluctuations can be observed at the early stage when 
the velocity fluctuations grow linearly in time. It becomes more significant for the small 
value of z Q /c, and can be estimated analytically from the ratio obtained from Eq. (THTT) as 

^| ^£^^o.4(^o/c) 2 , z n /c<l. (86) 

However, the values of ZqQo/c cannot be arbitrarily small, and are constrained by the un- 
derlying assumptions. To be consistent with the dipole approximation, the amplitude of the 
charged oscillator should be much smaller than its distance to the plate, and can be set to 
the order of 10 _2 zo ; for example. Thus, the corresponding velocity v of the charged oscillator 
is about v/c ~ 10~ 2 (z ^o/ c )- If can De argued that the charge's motion cannot be too slow 
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FIG. 2: Velocity fluctuations of the charge oscillator moving in the direction either parallel or 
normal to the plane boundary grow linearly at early times, and then saturate to a constant at 
late times. They start off at different rates, but approach the same saturated value h£l§j1m. The 
values of parameters are chosen as zq = zqQo/c = 1.0 and f c = r c Q,o/c = 5 x 10~ 4 , and e = 0.1; 
thus 7" 1 = 0.363 x 10 4 and j^ 1 = 0.931 x 10 4 . 

since it may give rise to a large position uncertainty to jeopardize the assumption of a point- 
like particle in the stochastic approach. It is of interest to take the electron as an example. 
Because the de Broglie wavelength, characterizing the size of the electron wavefunction, is 
^dB ~ h/(m e v) ~ 10~ 12 (i>/c) _1 m, when the velocity is overly small, the spatial extension of 
the electron wavefunction can be the same as or even larger than the distance to the plate 
Zq. Then it contradicts to the earlier assumptions. Nevertheless, the velocity of the electron 
may still be chosen as small as v jc ~ 10~ 4 , consistent with AdB ^ z ~ fim, ending up with 
z Q /c ~ 10~ 2 . Thus, the oscillation frequency for such a slow motion is reduced to the 
value of Qq ~ 10 12 s~ 1 . It will lead to the relaxation constant F ~ (e 2 f2 /m e c)fio ~ 10 s -1 
from Eq. (IHTT) . In addition, the dimensionless cutoff is then given by e = AdB^o/ c ~ 10~ 4 . 
Thus, even for an electron in rather slow motion, appreciable anisotropy Au 2 /At> 2 ~ 10~ 4 
can be found during the time regime 10~ 4 s < t < 0.1 s in which the cutoff dependence effect 
can be safely ignored while the saturation has not been reached yet in both directions. 
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At later times, the velocity fluctuations of the charged oscillator near the boundary reach 
saturation, independent of the orientation of the boundary, as shown in Fig. [2j The saturated 
value is approximately given by the Avf(oo) ~ M7 /2m as consequence of its nonuniform 
motion. The contribution of the velocity fluctuations, resulting from the imposition of the 
boundary conditions on electromagnetic fields, are at most the same order of magnitude 
as the cutoff-dependent terms. Both of them can be argued to be ignored at late times. 
Finally, the effective temperature corresponding to saturated velocity fluctuations for such 
Brownian motion can be estimated by 



where ks is the Boltzmann constant. 

To conclude, we study the influence due to quantized electromagnetic fields on the mo- 
tion of a nonrelativistic charged particle near the conducting plate. The nonlinear, non- 
Markovian Langevin equation of the particle is derived with the method of Feynman- Vernon 
influence functional, and it incorporates both dissipation backreaction on the charge in the 
form of the retarded self-force as well as the stochastic noise, manifested from vacuum 
fluctuations of quantized electromagnetic fields. The dipole approximation, an appropri- 
ate approximation for a non-relativistic motion, is implemented to find the solution to the 
Langevin equation. We consider that the charged particle undergoes a small-amplitude 
oscillation in the direction either parallel or perpendicular to the boundary plane. The 
noise-averaged trajectory of this charged oscillator is governed by the classical dynamics. 
Furthermore, the evolution of the kernel Ki(t), obtained by solving the Langevin equation, 
is found to be dominated by the narrow resonance in the weak coupling limit where the ratio 
of the decay width Tj over the oscillation frequency Qi is much less than unity, Ti/Qi <C 1. 
Thus, the memory effects or the the effects of the time-delay on backreaction terms due to 
the presence of the boundary are just to modify the quantities such as fij, Tj and so on. 
Then, velocity fluctuations of the charged oscillator driven by stochastic forces are found to 
grow linearly with time in the early stage of the evolution at the rate given by the relaxation 
constant. It turns out to be smaller in the parallel case than in the perpendicular one with a 
similar configuration, and reveals strong anisotropic behavior. They are then asymptotically 
saturated as a result of the fluctuation-dissipation relation at rather different relaxation time 
scales. However, we find that the same saturated value is obtained for motion in both direc- 




(87) 
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tions at late times, resulting from delicate balancing effects between dissipation backreaction 
and accompanying fluctuations, and thus the value is mainly determined by its oscillatory 
motion. So, at late times the effects from boundary-modified vacuum fluctuations on the 
velocity dispersion of the charged particle can be hardly seen. 

The dipole approximation amounts to linearizing the Langevin equation obtained above. 
However, beyond the dipole approximation, one expects to introduce the additional drift 
effects on the dynamics of the particle from the trajectory-dependent Green's functions. 
In particular, it may give rise to the noise-induced-drift forces owing to the correlations 
of stochastic forces. This effect has been considered in the context of the fast moving 



particle 



2l|, and may have observational consequences. 
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